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Abstract This contribution investigates situations in pedestrian dy- 
namics, where trying to walk the shortest path leads to largely differ- 
ent results than trying to walk the quickest path. A heuristic one-shot 
method to model the influence of the will to walk the quickest path is 
introduced. 

1 Introduction 

In 1952 John Wardrop raised his first principle that "The journey times 
in all routes actually used are equal and less than those which would be 
experienced by a single vehicle on any unused route." [1] which since has 
become a foundation of transport planning. But while all kinds of static 
and dynamic assignment methods [2-6] have been invented to calculate 
the user equilibrium of least journey time on road networks, in models 
of pedestrian dynamics [7] journey time mostly matters only indirectly. 
Typically the influence of following the shortest path is modeled explic- 
itly and the influence of the quickest path only indirectly with costs, 
reduced probabilities, or repelling forces, if an agent comes too close to 
another one, a group of pedestrians or an obstacle. These implicit meth- 
ods may have some effect for pushing agents more toward the quickest 
path. Additionally there are reasons, why journey time might matter less 
to pedestrians than for vehicle drivers, as attractiveness or safety of a 
route or just having plenty of time wanting to avoid the extra physical 
effort needed to walk on the quickest path may play a (more important) 
role. But considering the vast discussion of journey times in vehicular 
traffic, it appears that a discussion of the explicit influence of the will 
to minimize the expected remaining journey time is underrepresented 
in 2d models of pedestrian dynamics. Nevertheless, there is discussion 
on considering travel time and user equilibrium for network [8, 9] and 
macroscopic [10, 11] models of pedestrian flow. 

For the simulation of pedestrian dynamics things get more compli- 
cated, as basically there are two methods to consider journey time: either 



fully two dimensional without space semantic, respectively with full se- 
mantic symmetry between all accessible coordinates (this is discussed in 
[12]) or with a semantic decomposition and assignment of coordinates to 
special functional and semantic structures of the geometry i.e. the calcu- 
lation of a meaningful graph, which reduces the details of reality for the 
planning process. One may assume that when planning a path humans 
probably use something in between those alternatives: a two-dimensional 
representation of the space directly visible to them and a graph represen- 
tation of anything that's out of view. 

Generally, if the journey time itself is meant to be a determining factor 
for agents' behavior, there are only two methods of solution: analytically 
or iterative. There can not be a one-shot simulation considering journey 
time as input, as only at the end of the simulation the journey time is 
known. An analytic solution obviously is out of reach for any real-life 
application. An iterative solution is in principle possible, but in fact often 
impossible, if one considers typical computation times of simulations of 
pedestrian and evacuation dynamics. This leaves one-shot simulations 
using not journey time, but heuristics, and methods using simplifying 
assumptions like static assignment. 

It is obvious that the assumptions made for static assignment do not 
hold for evacuations. During evacuations there is not a constant, but a 
strongly time-dependent, short but strong demand, almost delta-function- 
like, widened by varying awareness times. During the egress time period 
of one individual the best decision - in terms of journey time - where 
to turn to at some decision spot, might change a number of times. The 
situation may be compared to the one of rush hours in vehicular traffic, 
and for these dynamic assignment methods were invented at first. 

Following this line of hand-waving argumentation, a heuristic seman- 
ticless method to model effects of the will to move on the path of least 
expected remaining journey time has been introduced [13] and will be 
sketched and investigated here. With reliance on a heuristic rather than 
estimated remaining travel time itself, in a strict sense the initial demand 
to consider the estimated travel time for walking behavior is given up. 
There are two arguments that may justify this: "estimation" in reality 
brings in a somewhat fuzzy element, such that real pedestrians also do 
not rely on a precise value of remaining travel time. It is not guaranteed, 
but it might even be that the heuristic discussed in this contribution might 
come closer to reality than a costly calculated remaining travel time. The 
second argument is that the results justify that it's more realistic to use 



the method than not to use it and if it's computationally faster than any 
method yielding the same benefit. 

2 Description of the Method 

Let's take a look at a group of pedestrians forming a jam at some bottle- 
neck, where alternative paths to detour around this group are still avail- 
able for other pedestrians approaching from behind. Assume the extra 
effort of walking the detour is reasonably small, such that a considerable 
ratio of real pedestrians would choose to walk it. How can it be modeled 
efficiently that agents of a simulation reproduce this behavior? 

In principle potentials [14-16] should do the task by just taking jams 
into accounts like obstacles. But problems arise with details: 

1. Some methods for potential calculation are not fast enough execution 
each time step. 

2. One has to deal with the case, when no detour is available, i.e. the 
jam would block the potential from spreading to some regions. 

3. The potential calculation method must be suited to consider the size 
of a jam quantitatively. A single agent must have only a negligible 
effect, a large jam should have a wide spreading effect. 

All this can be achieved using a flood fill, where a value of 1 is added, 
if a cell is unoccupied and some larger value s ac id, if it is occupied by 
an agent. But the two simple flood fill methods result in Manhattan or 
Chebyshev but not the desired Euclidean metric. In this sense these two 
methods produce large errors and unwanted artifacts in the movement. 
However, the errors can be reduced for once by combining the two meth- 
ods and second by using as heuristic not the time dependent potential, 
but its difference to the potential calculated on the empty (unoccupied) 
geometry. 

For two mutually visible coordinates (xo,yo) and (x±,y±), Manhattan 
distance is 

d M = \ Xl - x \ + \ Vl - y \ = \Ax\ + \Ay\ (1) 
and Chebyshev distance is 

<f = max.{\Ax\,\Ay\) (2) 

From this the "minimum distance" follows: 



d m = d M -d c = mm{\ Ax\,\Ay\) 



(3) 



and therefore the Euclidean distance is 



d E = ^J{d M ) 2 + {d m ) 2 = yJ\Ax\ 2 + \Ay\ 2 (4) 

Combining the whole potentials cell by cell in this manner - even if cells 
are not mutually visible - will now be called "method VI" . A discussion 
of the resulting errors of this method compared to Euclidean distance is 
given in [16]. 

Let's call the potential calculated on the empty (no agents) geometry 
Sy 1 and the potential after a certain time step Svi(t). Then the heuristics 
influencing the motion of pedestrians is Sd yn (t) = Syi(t) — S V1 . 

This has been used and coupled to the F.A.S.T. model [17-21] as 
an additional partial probability Pdyn = c 

k sdynS dy n(t) ^ multiplied to the 
original probability and normalized accordingly. 

3 Example of the Method's Effect 




Figure 1. Agents have to walk from the red to the green area, thereby 
choosing between one of the two corridors, of which one is 32 m longer. 

Figure 1 shows the scenario that is used in this contribution to in- 
vestigate the effect of the method. 4000 agents have to walk from the 
red to the green area. They have to choose one of the two corridors. The 
longer corridor is 32 m longer than the short one, which at the median 
maximum speed of the agents of v m = 1.6m/s equals to an additional 
20 seconds. Using the shorter corridor, the shortest total distance for an 
agent to walk is about 325 m, the longest about 418 m. This leads to a 
relative additional distance of 7.6% to 9.8%. 



Figure 2. Static potential with contrast and brightness adjusted to show 
the point inside the corridors, where the destination is equally far away, 
no matter, if an agent walks to the left or right. 

Figure 2 shows that there is no chance that without using a dynamic 
distance potential field an agent will move over the local maximum of 
distance toward the destination (i.e. local maximum of the static floor 
field) inside the longer corridor. Thus, in this case all agents will use the 
shorter corridor. 

For this scenario the method would cause, what is intended, if the load 
of the two corridors is shifted toward an equilibrium, implying a reduced 
total time for the process as well as a reduced average individual egress 
time of the agents. Figure 3 shows the values of these two observables in 
dependence of s a dd and the strong influence the method has, and figure 
5 directly shows the effect on route choice behavior. 

It is interesting to note that for s a dd = 2 and Isdyn = 1-0 not a 
single agent walks the longer corridor, still the two measured times drop 
significantly by about 12%. This is because the agents moved around 
the corners inside the corridor more efficiently. For higher values of s a dd 
more and more agents walk the longer corridor. This means that two time 
reducing effects are active. 

The dependence on parameter ksdyn surprisingly showed a more com- 
plicated behavior, which is shown in figure 4. Between kgdyn = 1-0 and 
ksdyn = 2.0 the evacuation times seem to stabilize and settle, just to start 
decreasing again at slightly higher values, pass a minimum and increase 
again and finally stabilize for very large values of ksdyn- 

A result as in figures 3 and 4 calls for a combined investigation of the 
two evacuation times' dependence on ksdyn and s a dd- The result of such 
an investigation for the average individual egress time is shown in figure 6. 
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Figure 3. Total times and average individual egress times in dependence 
of parameter s a dd- For each value of s a( id 100 simulations were carried out 
and the average was calculated. With increasing s a dd the total time and 
the average individual egress time slowly keep increasing, ksdyn has been 
set to ksdyn = 1-0 for these simulations. 
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Figure 4. Total times and average individual egress times in dependence 
of parameter k$dyn- Each value of kgdyn is an average of 100 simulations. 
For larger values of ksdyn the total time approaches 1020 seconds. s a dd 
has been set to s a dd = 10 for these simulations. 



Figure 5. Situation after 600 seconds without dynamic distance potential 
field (left) and with ksdyn = 1-0 and s a dd = 10 (right). 

The result for the total times is very similar, but the standard deviations 
are larger, as the arrival time of the last few agents varies increasingly 
with increasing product ksdyn " s add (the last few agents hesitate at the 
point of aequi-distance). 

In this two-dimensional parameter space the results look a bit complex 
as well. There are two local minima. (For larger values of the parameters, 
which cannot be displayed without obscuring the interesting region, the 
behavior qualitatively is the same as in figures 3 and 4.) 

At latest at this point one is interested in the number of agents walking 
through the longer corridor and the dependence of this number on the two 
parameters controlling the method. These numbers are shown in figure 7. 
With a maximum (at s a dd = 25 and ksd yn = 3.0) of 1543.8 agents using 
the longer corridor, the method achieves to do a considerable step toward 
the user equilibrium (which might be above 1900 agents). 

The second observation is that this plot just as well has two extrema 
(maxima in this case). This implies that the two maxima in figure 6 are 
not a result of the double influence of the method on local operational 
behavior (walking around a jam at a corner) and route choice behavior. 

On the other hand there's only a limited anti-correlation between the 
number of agents walking the longer corridor and the average individual 
egress time (see figure 8). This is a hint that increasing ksdyn a t some 
point leads to artifacts in the motion of the agents, which on average 
reduce the progress speed and therefore more agents need to walk the 
longer corridor to give the same average individual egress time as at 
lower ksdyn- I n table 1 pairs of results are compared, where either the 
numbers of agents walking the longer corridor are very similar and the 
average individual egress time is larger, where ksdyn 1S larger or vice 
versa the average individual egress times are similar and the number of 
agents in the longer corridor is smaller for smaller ksdyn- Interestingly for 
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Figure 6. Average individual egress times in dependence of k$dyn an d 

Sadd- 




Figure 7. Average number of agents walking the longer corridor in de- 
pendence of ksdyn and s a dd- Please note that this plot is rotated by 180 
degree compared to figure 6. 



identical kgdyn an d two different s a dd, there is no such effect, i.e. either 
both measured values are very similar or none. 



Correlation between # of agents in the longer corridor and the average individual egress time 

(averages within von Neumann neighborhood) 
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Figure 8. This figure shows the correlation between average individual 
egress time and load in the longer corridor of a data point and its nearest 
neighbor. Including the next to nearest neighbors (Moore neighborhood) 
does not change the result qualitatively, but makes it even more extreme 
quantitatively (the peak at ksdyn = 1-5 and s a( id = 25 is actually positive). 



4 Summary 



A non-iterational method was introduced that intends to shift the op- 
erational and tactic route choice behavior of agents in a simulation of 
pedestrian dynamics from a strong focus on the shortest path toward 
user equilibrium and a consideration of the quickest path. The investiga- 
tion with a proto-typic geometry gave promising results and showed that 
in general the idea works as originally thought off. 

However, in detail the results were more complicated than expected as 
dependencies on parameters showed up to be non-monotonic, yielding up 
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Table 1 . Comparison of similar results for different combinations of ksdyn 
and s a dd- "time" is the average individual egress time and "load" the 
number of agents passing the longer corridor. 



to three local extrema. The cause of this behavior could not be clarified 
entirely. 

A technical result is that values for parameter ksdyn that make agents' 
behavior look most natural do not realize the full potential of the method 
to come close to the user equilibrium. In turn this means that with param- 
eters lead to a route choice as close as possible to the user equilibrium, the 
operational behavior includes unrealistic artifacts from the construction 
of the dynamic distance potential field. 

It's in the nature of the method that the last agent leaving the longer 
corridor will always do this earlier than the last agent leaving the shorter 
corridor. But having these two agents doing this almost simultaneously 
is one criterion for Wardorp's first principle being fulfilled. I.e. with this 
method the last agent passing the shorter corridor would always have ar- 
rived earlier, if he had joined the last agent passing the longer corridor. 
Like already stated above: Wardrop's first principle is one thing, the ques- 
tion, how real pedestrians would behave, is different one. But if one takes 
the theoretical point of view, eager to have a simulation that leads to a 
result, where Wardrop's first principle is fulfilled, the method proposed 
in this contribution is a good starting point for an iterated dynamic as- 
signment search for an equilibrium. It will - in a computationally efficient 
way - reduce the number of iterations to find the equilibrium. 
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